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Background of Invention 

Field of the Invention 

[0001] The invention relates generally to well logging using nuclear magnetic resonance 

(NMR) instruments. More specifically, the present invention relates to methods and 
apparatus for NMR well logging based on carbon-proton J-coupling. 

Background Art 

[0002] Oil and gas exploration and production are very expensive operations. Any 

knowledge about the formations that can help reduce the unnecessary waste of 
resources in well drilling will be invaluable. Therefore, the oil and gas industry have 
developed . various tools capable of determining and predicting earth formation 
properties. Among different types of tools, nuclear magnetic resonance (NMR) 
instruments have proven to be invaluable. NMR instruments can be used to determine 
formation properties, such as the fractional volume of pore space and the fractional 
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volume of mobile fluid filling the pore space. General background of NMR well 
logging is described in U.S. Patent No. 6, 140,817. 

[0003] Nuclear magnetic resonance is a phenomenon occurring in a selected group of 

nuclei having magnetic nuclear moments, i.e., non-zero spin quantum numbers. When 
these nuclei are placed in a magnetic field (B 0 , "Zeeman field"), they each precess 
around the axis of the B 0 field with a specific frequency, the Larmor frequency (a) 0 ), 
which is a characteristic property of each nuclear species (gyromagnetic ratio, y) and . 
depends on the magnetic field strength (B Q ) effective at the location of the nucleus, i.e., 
(Do - Y B 0 . 

[0004] Proton is the major nucleus of investigation in well logging NMR applications 

because of its good NMR sensitivity and its high abundance in water and hydrocarbons. 
Furthermore, due to downhole limitations, the current well logging tools only measure 
fi, T2 relaxation times and diffusion effects 

[0005] In other fields of NMR applications such as chemistry, biology, and petroleum 

fluid analysis, proton and carbon chemical shift and J-coupling spectroscopic 
techniques are routinely used to determine molecular structures. Chemical shift is the 
term given to describe the screening effect of the electrons to the magnetic field that a 
nucleus experiences. Different chemical groups, such as CH2 and CH3, have different, 
magnitude of screening effects, and, therefore, they appear as separate peaks in the 
proton chemical shift spectrum. The separation in frequency of different peaks is 
proportional to the static magnetic field strength, i.e., magnetic field dependent. J- 
coupling, also known as spin-spin or scalar coupling, originates from spin interaction 
between nuclei through bonding electrons. See, E. L. Hahn, and D. E. Maxwell, Spin 
echo measurements of nuclear spin coupling in mokcules y Physical Review 88, 1 070- 
1084 (1952). J-coupling experiment is seldom performed by itself Instead, J couplings 
are always measured together with chemical shifts through one-dimensional or multi- 
dimensional spectroscopic techniques. 

[0006] As noted above, chemical shift is magnetic field dependent. The homogeneity of 

a static magnetic field has to be within a few parts per million (ppm) to perform 
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chemical shift spectroscopy. This level of homogeneity is hard to realize in the 
wellbore using existing technologies. In contrast, J-coupling constants are independent 
of static magnetic field strengths and temperatures. This makes it possible to perform J- 
coupling experiment without chemical shift spectroscopy in inhomogeneous static and 
radio-frequency (RF) magnetic fields. 

[0007] U.S. Patent No. 6,111,409 issued to Edwards discloses methods for performing 

chemical shift spectroscopy in the wellbore. Because the static magnetic field 
homogeneity is difficult to achieve in the formation, the methods of Edwards involve 
withdrawing fluids into a formation tester before NMR measurements. Even in the 
formation tester, the stringent homogeneity required for conventional chemical shift 
measurements is not an easy task. A permanent magnet and shim coils will be used to 
generate the static magnetic field. First of all, the space in the tool that can 
accommodate the magnet and shim coils is very limited in shape and size. To achieve a 
1 ppm homogeneity over a reasonable volume is a daunting task for magnet design and 
manufacturing. If the homogeneous volume is too small, the small sample may not be a 
good representation of the fluid being investigated. Secondly, temperature change can 
affect the strength and homogeneity of the magnetic field. 

[0008] U.S. Patent No. 6,346,813 issued to Kleinberg discloses a variety of NMR 

measurements for characterizing fluid samples withdrawn from subsurface formations. 
This patent is assigned to the assignee of the present invention and is hereby 
incorporated by reference. Proton and carbon NMR measurements are among the 
proposed techniques. 

[0009] U.S. Patent Application Serial No. 10/064,529 filed by Speier on July 24, 2002, 

discloses methods for J-spectroscopy experiments using spin-echo difference 
techniques. This application is assigned to the assignee of the present invention and is 
hereby incorporated by reference. This approach provides convenient methods to 
obtain J couplings in the wellbore. Although this approach is less sensitive to magnetic 
field inhomogeneity, it still suffers from field inhomogeneity to some extent. The 
inhomogeneous) B 0 (static) field makes the effective Bi (radio-frequency) field 
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inhomogeneous. As a result, the 71 pulses may not be accurate throughout the region of 
investigation. Inaccurate 71 pulses may diminish the J-modulation signals. 

[0010] Therefore, it is desirable to have NMR methods and apparatus for determining J 

coupling that are less affected by magnetic field inhomogeneity. 

Summary of Invention 

[0011] In one aspect, embodiments of the invention relate to methods for performing 

nuclear magnetic resonance measurements having J modulation information. A method 
for obtaining nuclear magnetic resonance measurements in accordance with one 
embodiment of the invention includes inducing a static magnetic field in a formation 
fluid sample; applying an oscillating magnetic field to the fluid sample according to a 
preparation pulse sequence that comprises a J-edit pulse sequence for developing J 
modulation; and. acquiring the nuclear magnetic resonance measurements using a 
detection sequence, wherein the detection sequence comprises at least one 180-degree 
pulse. The method may further include acquiring the nuclear magnetic resonance 
measurements a plurality of times each with a different value in a variable delay in the 
J-edit pulse sequence; and analyzing amplitudes of the plurality of nuclear magnetic 
resonance measurements as a function of the variable delay to provide J coupling 
information. 

[0012] Another aspect of the invention relates to methods for characterizing formation 

fluids. A method for characterizing formation fluids in accordance with one 
embodiment of the invention includes providing a nuclear magnetic resonance 
instrument in a borehole; inducing a static magnetic field in a region of interest; 
applying an oscillating magnetic field to the region of interest according to a 
preparation pulse sequence that comprises a J-edit pulse sequence for developing J 
modulation; and acquiring nuclear magnetic resonance measurements using a detection 
sequence, wherein the detection sequence comprises at least one 180-degree pulse. 

[0013] Another aspect of the invention relates to nuclear magnetic resonance instruments 

for well logging. A nuclear magnetic resonance instrument in accordance with one 
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embodiment of the invention includes a housing adapted to move in a wellbore; a 
magnet disposed in the housing adapted to induce a static magnetic field in a zone of 
interest; an antenna assembly disposed in the housing, the antenna assembly adapted to 
induce an oscillating magnetic field in the zone of interest and to receive nuclear 
magnetic resonance signals; and an electronic module including a memory to store 
instructions for performing a J-edit pulse sequence 

[0014] Other aspects of the invention will become apparent from the following 

. description, the drawings, and the claims. 

Brief Description of Drawings 

[0015] FIG. 1 shows a schematic diagram of an NMR logging tool in a well bore suitable 

for practicing methods of the invention. 

[0016] FIG. 2 shows a schematic diagram showing a pulse sequence for J editing 

according to one embodiment of the invention. 

[0017] FIG. 3 shows a graph of signal amplitudes as a function of the delay time for J 

modulation according to one embodiment of the invention. 

[0018] FIG. 4 shows a schematic diagram showing a pulse sequence for J editing 

according to one embodiment of the invention. 

[0019] FIG. 5 shows a schematic diagram showing a pulse sequence for J editing 

including gated decoupling according to one embodiment of the invention. 

[0020] FIG. 6 shows a schematic diagram showing a pulse sequence for J editing 

including gated decoupling according to one embodiment of the invention. 

[0021] FIG. 7 shows a schematic diagram showing a pulse sequence for J editing with 

proton detection according to one embodiment of the invention. 

[0022] FIG. 8 shows a schematic diagram showing a pulse sequence for J editing with a 

signal enhancement pulse sequence according to one embodiment of the invention. 
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[0023] FIG. 9 shows a schematic diagram showing a pulse sequence for J editing with a 

signal enhancement pulse sequence and gated decoupling according to one embodiment 
of the invention. 

[0024] FIG. 10 shows a schematic diagram showing a pulse sequence for J editing with 

magnetization transfer from proton to carbon according to one embodiment of the 
invention. 

[0025] FIG. 1 1 shows a schematic diagram showing a pulse sequence for J editing with 

magnetization transfer from carbon to proton according to one embodiment of the 
invention. 

[0026] FIG. 12 shows a schematic diagram showing a pulse sequence for J editing with 

magnetization transfer from proton to carbon and then back to proton according to one 
embodiment of the invention. 

[0027] FIG. 13 shows a flow chart of a method for deriving J coupling and for 

hydrocarbon typing according to one embodiment of the invention. 

[0028] FIG. 14 shows a flow chart of a method for calculating oil- water ratios and for 

fluid typing according to one embodiment of the invention. 

Detailed Description 

[0029] Embodiments of the invention relate to apparatus and methods for NMR well 

logging based on J-couplings (e.g., ^-^C). Embodiments of the invention are based on 
measuring the effects of J-coupling (e.g., ! H- 13 C) on spin echo amplitudes. Methods in 
accordance with the invention are less sensitive to magnetic field inhomogeneity. 
Modulation of spin echo amplitudes by the J-coupling effects is referred to as "J- 
editing" in this description. J-editing can be used to determine detailed compositions of 
an oil sample, for example, the abundances of various carbon groups. These carbon 
groups may include quaternary C, aromatic CH, aliphatic CH, CH 2 , CH 3 , and CH 4 . 

[0030] Methods of the invention may be practiced with an NMR instrument similar to 

those known in the art. The NMR instrument may be in the laboratory or be part of a 



5 



PATENT APPLICATION 
ATTORNEY DOCKET NO. 20.2862 (R&O 09428.106001) 

wireline logging tool, a measurement-while-drilling (MWD) tool, or a logging-while- 
drilling (LWD) tool. In addition, the NMR instrument may be part of a formation tester 
(fluid sampling tool) used to analyze the fluids withdrawn from the formations. 

[0031] FIG. 1 illustrates a schematic of an NMR well logging tool 30 placed in a 

wellbore 32 that penetrates a formation 31. The NMR logging tool 30 is suspended in 
well bore 32 on an armored cable 33, the length of which substantially determines the 
relative depth of the tool 30. The cable length is controlled by suitable means at the 
surface such as a drum and winch 8. Surface equipment 7 may comprise a processor 
which communicates with the downhole equipment. Although a wireline tool is shown 
in FIG. 1, variations of NMR logging tools may include MWD or LWD tools. The tool 
30 may include one or more permanent magnets, represented as block 34, and one or 
more antennae (e.g., an antenna assembly), represented as block 35, that can induce RF 
(or oscillating) magnetic field and function as receivers. In addition, the tool 30 may 
include an electronic module 36, which may include a central processor 36a, a memory 
36b, and other components (not shown) for controlling the pulse sequences and the 
acquisition. The memory 36b may store instructions for performing the J-edit pulse 
sequence in accordance with embodiments of the invention. 

[0032 J J-editing is based on the presence of J couplings between different nuclei, such as 

between different protons (homonu clear coupling) and between carbons and protons 
(heteronuclear coupling). See e.g., L. Muller, A. Kumar, and R. R. Ernst, Two- 
dimensional carbon-13 NMR spectroscopy, J. Chem. Phys. 63, 5490-5491 (1975). 
While methods of the invention may be applied to both homo- and hetero-nuclear 
couplings, for clarity, the following description uses only hetero-nuclear coupling 
experiments as illustrations. In addition, the following description uses carbon-proton 
coupling as an example. However, one of ordinary skill in the art would appreciate that 
methods of the invention may also be applied to other hdmonuclear or heteronuclear 
couplings. 

[0033] A method in accordance with the invention generally comprises two parts. In the 

first part, the precession of observed spins (e.g., carbon) due to J-coupling to coupled 
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spins (e.g., proton) is produced, resulting in an NMR signal of observed spins (e.g., 
carbon) that is a function of the J-coupling and some parameters of the pulse sequence. 
The pulse sequences used in the first part for developing J coupling effects are referred 
to as "J-edit pulse sequences" in this description. In the second part, the J-modulated 
observed spin (e.g., carbon) signal is detected, preferably using a pulse sequence (e.g., a 
% pulse train as in a Carr-Purcell-Meiboom-Gill, CPMG, sequence) that maximizes 
signal-to-noise ratios (SNR), The pulse sequences used in the second part are referred 
to as "detection sequences." One of ordinary skill in the art would appreciate that 
various detection sequences, including CPMG-like sequences, may be used without 
departing from the scope of the invention. Thus, a complete pulse sequence in 
accordance with some embodiments of the invention comprises a J-edit pulse sequence 
followed by a detection sequence. In some embodiments of the invention to be 
described later, a complete pulse sequence may further include a signal enhancement 
pulse sequence (e.g., nuclear Overhauser effect (NOE) or magnetization transfer pulse 
sequences) before or after the J-edit pulse sequence. The pulse sequences before the 
detection pulse sequence is generally referred to as a "preparation pulse sequence " 
Thus, a preparation pulse sequence may include only a J-edit pulse sequence or 
additionally include a signal enhancement pulse sequence. 

[0034] In J-editing experiments, radio frequency (RF) pulses at both carbon and proton 

resonance frequencies are generated and transmitted to the probe. The probe can use 
either one or two RF coils to induce and detect proton and carbon NMR signals. If two 
coils are used, one coil, is tuned to the proton resonance frequency and the other to the 
carbon frequency. If one coil is used, the coil is double-tuned to both proton and carbon 
resonance frequencies. For one example, see F. D. Doty, R. R. Inners, and P. D. Ellis, A 
multinuclear double-tuned probe for applications with solids or liquids utilizing lumped 
tuning elements, J. Magn. Reson. 43, 399-416 (1981). Those of ordinary skill in the art 
would appreciate that other techniques may be used to induce and detect signals. An 
antenna assembly is used to refer to such coils, whether a single coil or two coils. 

[0035] J-editing methods in accordance with the invention may be categorized into two 

basic groups - carbon detection methods and proton detection methods. In carbon 
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detection methods, 13 C spins are the source of magnetization and the subject of 
detection (i.e., the observed spins), and the proton spins (i.e., the coupled spins) are 
flipped or decoupled to impose a J-moduIation on the 13 C spin echo amplitude. This is 
conventionally denoted as 13 C{ 1 H}, in which the nucleus before the bracket is the 
observed nucleus and the nucleus in the bracket is the coupled/decoupled nucleus. In 
proton detection methods, proton spins are the source of magnetization and the subject 
of detection (i.e., the observed spins), and the 13 C spins (i.e., the coupled spins) are 
flipped or decoupled to modulate the proton spin echo amplitude through J-coupling. 
This is conventionally denoted as^ ^{^C}. The signal contribution from protons 
attached to 12 C is not modulated by the j-coupling and, therefore, may be removed by 
subtracting the signals of two scans, one with J-modulation and the other without. 

[0036] The NMR sensitivity of 13 C detection is low compared to that of proton 

detection because of the low gyromagnetic ratio (y) and the low natural abundance 
(about 1.1%) of the 13 C species. This makes it desirable to enhance the sensitivity of 
the carbon detection J-editing experiments. A variety of signal enhancement techniques 
may be used for this purpose. The signal enhancement techniques that may be used 
with embodiments of the invention, for example, may include nuclear Overhauser 
enhancement (NOE) techniques and magnetization transfer techniques. See G. A.- 
Morris, and R. Freeman, Enhancement of nuclear magnetic resonance signals by 
polarization transfer, J Am. Chem. Soc. 101:3, 760-762 (1979); S. L. Hartmann and E. 
L. Hahn, Nuclear double resonance in the rotating frame, Phys. Rev. 128, 2042-2053 
(1962). 

J-Editing with Carbon Detection 

[0037] Carbon detection J-editing methods have two subgroups: proton flip methods and 

gated decoupling methods. These methods will be described in more detail in the 
following sections. 

[0038] FIG. 2 shows a method according to one embodiment of the invention. This 

method uses a proton flip pulse in the decoupler channel to develop J modulation in the 
observation (e.g., 13 C) channel. The pulse sequence as shown in FIG. 2 may be 
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separated into two parts. The first part comprises a preparation pulse sequence that 
includes only a J-edit pulse sequence during the J-modulation period. The J-edit pulse 
sequence includes a typical spin-echo ((7i/2)-delay-(7t)-delay) pulse sequence in the 13 C 
channel and a k pulse in the ! H channel. The second part comprises a detection 
sequence that includes a train of n pulses, as in a typical Carr-Purcell-Meiboom-Gill 
(CPMG) pulse, for the measurement of the spin echoes. Such a n pulse train in the 
detection sequence is generally referred to as a CPMG-like sequence in this description. 

[0039] As shown in FIG.-2, the first echo spacing for carbon (t e i) is usually large, e.g., on 

the order of 1/J » 5-8ms, such that the J modulation would have sufficient time to 
develop. The later echoes (t e ) in the CPMG-like sequence preferably have short echo 
spacing (e.g. 200 )as) so .that a large number of echoes may be acquired to improve SNR 
before diffusion and T 2 decays reduce signal amplitudes. 

[0040] As shown in FIG 2, a single n pulse is applied in the proton channel at time x 

after the first n/2 carbon pulse. Because of the C-H J coupling, after the initial n/2 
carbon pulse, the carbon magnetization splits into two vectors preeessing around the 
static magnetic field B 0 at two different angular frequencies. The changing phase angle 
(separation) between the two vectors causes the total magnetization amplitude to change 
as a cosine function of the phase angle. In other words, the amplitude of the carbon 
magnetization is J-modulated. During the x period, the phase angle (separation) 
between the two vectors gradually increases, until the proton n pulse exchanges the 
angular frequencies of the two vectors. As a result of the proton n pulse, which 
exchanges the spin states of the proton spins, the phase angle between the two vectors 
decreases until the carbon n pulse is applied. The carbon n pulse exchanges the spin 
states of the carbon spins, and hence the angular frequencies of the two vectors are 
exchanged again. Thus, the phase angle (i.e., phase separation between the two vectors) 
increases again. From this description, it is clear that J modulation only occurs during a 
period corresponding to 2x within the first echo spacing. In other words, the extent of J 
modulation depends on the delay time x. This J modulation manifests itself as a J- 
modulation factor that modulates the first echo amplitude. 
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[0041] After the first echo, only carbon n pulses are applied to the carbon channel to 

generate an echo train, but no proton pulses are applied to the proton channel; As a 
result, the carbon % pulses, which are applied in the middle of. the echo spacings, 
suppress further J modulation. As a result, J modulation occurs only during the first 
echo spacing, i.e., the J-modulation factor is frozen after the first echo and all 
subsequent echoes have the same J modulation as the first echo. Because all echoes 
have the same J-modulation factor, it is possible to sum up all echoes to improve the 
signal-to-noise ratio (SNR) without losing J-modulation information. When all echoes 
are summed up, the T2 and diffusion decay effects in the second and subsequent echoes 
are ignored. This has little effect on the derivation of the J-modulation factor. 
Alternatively, the J-modulation factor may be found by a fitting procedure where T 2 and 
diffusion decay effects in the second and subsequent echoes are taken into account. 

[0042] As noted above, J modulation is frozen after the first echo. Therefore, the J- 

modulation factor may be derived from the first echo amplitude S(x). Without 
considering the effects of static and RF field inhomogeneities, S(x) is given by 

S(x) = S c 0 Y J A l (cos2%J ) T) ni Eq. (1) 

i 

where S c 0 is the carbon spin echo amplitude without the J modulation, Ai is the relative 
abundance of the i ih carbon group, ] l the proton-carbon J-coupling constant for the i th 
carbon group, and nj is the number of protons in the i th carbon group. 

[0043] Eq. (1) shows that with a fixed x, the amplitude of the first echo (or the sum of all 

echoes) is related to cosine functions of the coupling constant (Jj) to the n r th power, 
where n, is the number of protons attached to the carbon. The number of protons (nj) 
attached to a carbon may have five different values: 0 for C (quaternary carbon), 1 for 
CH, 2 for CH 2 , 3 for CH 3 , and 4 for CH 4 . The number of protons (n* = 1-4) attached to 
a carbon has little effect on the C-H J-coupling constants. For example, C-H coupling 
constants for all aliphatic carbons (CH, CH 2 , CH3, and CH4) are around 125 Hz 

[0044] In Eq. (1), carbon groups are defined in such a way that different carbon groups 

either have different J-coupling constants (Jj) or different number of protons (nj). 
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Because J-cpupling constants do not vary much, it is convenient to group the J-coupling 
constants into a few values. For example, all the J-coupling constants in an oil sample 
may be approximately grouped into two values: 160 Hz for aromatic carbon groups and 
125 Hz for aliphatic carbon groups. By doing so, the oil sample may have six different 
carbon groups, which are listed in table 1. For typical oil samples, there is little 
aliphatic CH and most of the CH is in aromatics. Thus, the aliphatic CH group may be 
ignored to simplify the model, if desired. 

[0045] Table 1. An example of carbon group modeling 



Carbon group 


Name 


C-H coupling 


Number of 


index i 




constant Jj 


protons ni 


1 


Quaternary C 


None 


0 


2 


Aromatic CH 


160 Hz 


1 


3 


Aliphatic CH 


125 Hz 


1 


4 


CH 2 


125 Hz 


2 


5 


CH 3 


125 Hz 


3 


6 


CH4 


125 Hz 


4 



[0046] With a particular sample, the distributions of Ji and nj do not change, if a series of 

measurements are acquired with different delay times, the amplitudes of the first echoes 
in these measurements would depend on cosine functions of the delay time x. 
According to a preferred embodiment of the invention, multiple scans with different t 
values may be acquired with a fixed first echo spacing such that the T2 and diffusion 
effects are the same for all measurements. When these measurements are acquired with 
different delay times t, the amplitudes of these measurements would depend on the 
delay times t. The signal amplitudes as a function of x would form a set of linear 
equations, from which A;, the abundances of different carbon groups, may be computed. 

[0047] In solving these linear equations, the J-coupling constants may be obtained from 

the published values in the literature or from lab NMR spectroscopic experiments on oil 
samples. As noted above, the J coupling constants do not vary significantly among all 
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aliphatic carbon groups and these coupling constants can be approximated with a few 
values that depend on the hybridization of the carbon (e.g., the values shown in Table 
1). 

[0048] In the above described approach, a set of linear equations may be obtained with a 

plurality of different x values, which need not have any particular relationship to each 
other. In an alternative approach, a plurality of measurements may be acquired with 
uniformly incremented x values. With such a data set, Fourier transform may be used to 
produce a J spectrum, from which J; and A\ may be derived. 

[0049] FIG. 3 shows an example graph of echo amplitudes as a function of t delays using 

a light oil sample. This example includes 3 1 x values from 0-8 ms. The modulation of 
the echo amplitudes by the x delays is apparent from this graph. It is clear from this 
graph that a majority of the echo amplitudes is modulated by a J coupling that manifests 
itself as a periodicity of about 8 ms in this graph. This periodicity corresponds to a J 
coupling constant of about 125 Hz, indicating that the majority of signals are from 
aliphatic carbon groups. This is consistent with the composition of the sample. 

[0050] FIG. 2 illustrates only a basic scheme for a method in accordance with the present 

invention. One of ordinary skill in the art would appreciate that other techniques known 
in the art may be included to improve the measurements. For example, phase cycling 
may be included to reduce or suppress constant noises (e.g., DC offset of the 
instrument) or signal degradation due to imperfect pulse lengths (i.e., the n/2 or n 
pulse). Thus, the use of phase cycling is expressly within the scope of the present 
invention. Examples of the phase cycling schemes include CYCLOPS, which can be 
used to remove DC offset and unbalanced gains in quadrature detection. In addition, 
cycling of the first n/2 pulse can remove the extra carbon magnetization recovered due 
to T] relaxation during the long train of CPMG-like k pulses. 

[0051] Furthermore, one of ordinary skill in the art would appreciate that the 7c/2 and % 

pulses depicted in FIG. 2 may be substituted with composite pulses. Various composite 
k/2 and n pulses are known in the art. See M.H Levitt and R. Freeman, ./. Mag. Reson. 
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(33), 473 (1979). For example, a composite four-pulse sequence, 
P ^/j j^^j ^/j {x/^j 9 is commonly used as a substitute for a 7i/2 pulse in the X- 

axis direction, (n/2) Xi See M.H. Levitt, J. Mag. Reson. (48), 234 (1982). Similarly, a 
7i pulse, may be substituted by, for example, the following composite pulse: 

^ = (/^) : T ^ ese °P m P os i te rc/2 or pulses typically produce better results 

because they are less prone to signal degradation due to imperfect pulse lengths for the 
k/2 or n pulses. The pulse length imperfection is expected to be more severe in the 
downhole environment, where the temperature is expected to have a significant impact 
on the characteristics of the intricate circuitry in the pulse programmer and transmitter. 
Note that the use of composite pulses is expressly within the scope of the present 
invention 

[0052] As shown in FIG. 2, the proton n pulse is applied after a delay time x. This 

proton k pulse does not need to go past the first carbon n pulse (i.e., x < Ui/2). 
However, embodiments of the invention are not so limited. In some embodiments, the 
proton 7i pulse may coincide with the first carbon % pulse (see FIG. 4). In other 
embodiments, the proton n pulse may go past the first carbon n pulse. In this case, x is 
defined as the time delay between the proton n pulse and the first echo such that Eq. (1) 
still holds. 

[0053] FIG 4 shows a pulse sequence for another proton flip method in accordance with 

embodiments of the invention. In this pulse sequence, the proton n pulse and the 
carbon % pulse are applied simultaneously or with a fixed time shift. The first echo 
spacing and x are varied to produce different J modulations. In this case, the first echo 
spacing is defined as 2x and Eq. (1) still holds. With this pulse sequence, the T2 and 
diffusion effects will be different when x (hence the first spacing) is varied. Therefore, 
if a series of measurements are acquired with a plurality of x values, correction may 
need to be applied to remove/correct the T 2 and diffusion effects in these measurements 
before these data are analyzed according to methods described above (e.g., solving a set 
of linear equations or performing a Fourier transformation). 
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[0054] As shown in FIG. 2 and FIG. 4, a pair of carbon/proton n pulses are applied 

during the first echo spacing, but only carbon n pulses are applied during the 
subsequent echo spacings in the CPMG-like pulse train. With this pulse sequence, J 
modulation builds up during the first echo spacing only. One of ordinary skill in the art 
would appreciate that other modifications to these pulse sequences are possible without 
departing from the scope of this invention. For example, it is also possible to build up 
the J-modulation during the first few echo spacings. During each of these early echo 
spacings, which should be relatively long (e.g., 8 ms) for J modulation to develop, a pair 
of carbon and proton it pulses are applied, as shown for the first echo spacing in FIG. 2 
and FIG. 4. One of ordinary skill in the art would also appreciate that this technique of 
using multiple pairs of carbon/proton n pulses to build up J- modulation in the first few 
echo spacings may be applied to all proton-flip carbon-detection methods and carbon- 
flip proton-detection methods which will be described later. 

[0055] J modulation may be created with n pulses (flip pulses) or with gated decoupling 

in the coupled nuclei channel. Some embodiments of the invention use gated 
decoupling, instead of flip pulses. FIG. 5 shows one such pulse sequence, in which a 
proton decoupling replaces the proton n pulse in FIG. 2. The proton decoupling stops J- 
modulation buildup from time x to the first echo. Thus, J modulation develops only 
during t period using the pulse sequence of FIG, 5 as opposed to 2x period using the 
pulse sequence of FIG. 2. As a result, the first echo amplitude S(x) from the experiment 
shown in FIG. 5 is given by: 

S^^SSIA^coskJ^ Eq. (2) 

i 

[0056] In a related pulse sequence, the proton decoupling may be applied right after the 

carbon n/2 pulse and be removed at some point before the first echo. In this case, x is 
defined as the time delay between the end of the decoupling pulse and the first echo and 
Eq. (2) still holds. One of ordinary skill in the art would appreciate that proton 
decoupling as depicted in FIG. 5 may be accomplished by any composite pulses known 
in the art for decoupling or with a broad-band (e.g., noise modulation) decoupling pulse. 
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In addition, while FIG. 5 shows that C is the observed spin and H is the coupled spin 
(i.e., 13 C{ 1 H}) for clarity of illustration, these pulse sequences may also be applied to 
other hetero nuclear combinations. 

[0057] FIG. 6 shows another embodiment of a gated-decoupling method. As shown in 

FIG. 6, the proton decoupler covers either the first half or the second half of the first 
echo spacing In addition, the first echo spacing is varied to produce different J 
modulations, which are encoded in the different amplitudes of the recorded signals. 
This approach is similar to that shown in FIG. 4. A series of spectra may be recorded 
and the modulation of signal amplitudes as a function of the first echo spacing (or t) 
may be analyzed (e.g., with a set of linear equantions or with Fourier transformation) to 
provide the J coupling constants and Aj. 

J-Editing with Proton Detection 

[0058] The above described embodiments are based on carbon detection and proton 

flipping or decoupling, i.e., 13 C{ l H}: As noted above, carbon detection is much less 
sensitive compared with proton detection, due to the lower gyromagnetic ratio (y) and the 
lower abundance of the carbon spin. Therefore, it is sometimes desirable to perform 
these experiments with proton detection, i.e., ^{^C}. 

[0059] FIG 7 shows a proton detection pulse sequence similar to that shown in FIG. 2, 

except that the l3 C and l H channels are swapped. In this proton detection pulse 
sequence, proton spins are the source of magnetization and the subject of detection (i.e., 
the observed spins). Carbon spins are flipped to enable the buildup of the J-modulation. 
Because about 99% of the proton signal comes from protons attached to 12 C and has no 
J coupling, the signal contribution from protons attached to 12 C should be removed to 
facilitate analysis. . This may be accomplished by subtracting the signals of two 
measurements, one with the carbon % pulse (i.e., J modulation) and the other without. 
Neglecting the field inhomogeneity effects, the first echo amplitude S(t) after the 
subtraction is given by: 
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S(T) = ^H^XA,n,(l-cos2 n J,T) Eq. (3) 

i . 

where the factor 0.011 is the natural abundance of 13 C, s{| is a proton spin-echo 
amplitude without J modulation. The terms A i? j ; , and n; are as defined in Eq. (1). 

[0060] In carbon detection methods as illustrated in FIGs. 2 and 4-6, n 4 changes J- 

modulation frequency because each carbon is simultaneously influenced by nj protons 
directly attached to it. In contrast, with a proton detection method as illustrated in FIG. 
7, nj does not change J-modulation frequency because each proton is only influenced by 
one carbon directly attached to it. As a result, the proton-detection J-editing method 
cannot distinguish different carbon groups based on their differences in n^ However, 
this method can still differentiate different carbon groups if they have different J- 
coupling constants, e.g., aromatic C-H coupling » 160 Hz, and aliphatic C-H coupling « 
125 Hz. Thus, it is feasible to use this proton-detection J-editing method to distinguish 
aliphatics and aromatics based on different J's. 

[0061] From Eq. (3), it is apparent that the signal contribution from a certain carbon 

group reaches its maximum value when 

x = l/2J i . Eq. (4) 

It is also known that the J-coupling constants do not vary much among different carbon 
groups in oil samples, which contain mostly aliphatic hydrocarbons. Therefore, a 
/ properly selected (average) t can be used to ensure that signal contributions from all 

carbon groups reach levels close to their maximum values. 

[0062] One advantage of this approach is that the measured signal amplitude S(x) is a 

good indication of the hydrogen index of hydrocarbons even in the presence of water in 
the sample. Thus, the measured signal amplitude S(x) using the pulse sequence shown 
in FIG. 7 together with a regular proton measurement can be used to provide a 
convenient method for computing the oil/water ratio in the sample. Accordingly, this 
approach provides a good fluid typing method. 
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[0063] FIG. 7 illustrates only one exemplary pulse sequence for proton detection based 

on the pulse sequence shown in FIG. 2. One of ordinary skill in the art would 
appreciate that different proton detection J-editing pulse sequences may also be devised 
that are analogous to pulse sequences illustrated in FIGs. 4-6. Specifically, a pulse 
sequence for the proton-detection method with varied echo spacing can be obtained by 
swapping the *H and 13 C channels as illustrated in FIG. 4. Similarly, gated-decoupling 
proton-detection J-editing pulse sequences may be obtained by swapping the *H and l3 C 
channels illustrated in FIG. 5 and FIG. 6. 

Signal Enhancement Techniques 

Nuclear Overhauser Enhancement 

[0064] As noted above, with carbon detection the sensitivity is low. In some methods 

according to the invention, the carbon detection pulse sequences may be combined with 
signal enhancement pulses. Various techniques are known in the art for enhancing 
signal amplitudes from insensitive nuclei. These techniques, for example, include 
Nuclear Overhauser Effects (NOE) and polarization (magnetization) transfer. 

[0065] With NOE enhancement techniques, the signal amplitudes for the carbon 

detection J-editing methods may be enhanced by applying a proton broadband 
decoupler, before carbon experiments. Such NOE enhancement can produce up to three 
times signal enhancement because the gyromagnetic ratio of a proton is about four 
times that of a carbon. The broad-band decoupling may be accomplished with an RF 
transmitter with noise modulation. Alternatively, this may be accomplished with a 
series of pulses with very short spacings. The duration and phase of each pulse in the 
series may be varied to generate broadband decoupling. 

[0066] FIG. 8 shows the pulse sequence of FIG. 2 combined with a NOE pulse sequence, 

while FIG. 9 shows a gated-decoupling J-editing pulse sequence (as shown in FIG. 5 or 
FIG. 6) incorporated with an NOE pulse sequence. These pulse sequences are only for 
illustration. One of ordinary skill in the art would appreciate that NOE decoupler/pulse 
sequences may be applied to all carbon detection J-editing methods. As shown in FIGs. 
8 and 9, the signal enhancement pulse sequence (the NOE pulses) together with the J- 
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Edit pulse sequence are referred to as a preparation pulse sequence. A preparation pulse 
sequence as used herein refers to pulse sequences prior to the detection sequence. Thus, 
in some embodiments, the preparation pulse sequence may include only the J-Edit pulse 
sequence, as shown in FIG. 7, while in other embodiments, the preparation pulse 
sequence may include both the signal enhancement pulse sequence and the J-edit pulse 
sequence, as shown in FIGs. 8 and 9. 

[0067] Another approach to signal enhancement uses polarization transfer (or 

magnetization transfer). For example, hetero-nuclear magnetization transfer pulse 
sequences may be incorporated into the carbon detection J-editing methods to improve 
SNR. The direction of the magnetization transfer can be 'H-^C, l3 C-VH, and l H 
-V 3 C-Vh. The SNR of the detected signal for NMR experiments is given by 

SNRocy SOUfce ( Yobserve ) 3/2 Eq.(5) 

where y SO urce is the gyromagnetic ratio of the source spins, y 0 bserve the gyromagnetic 
ratio of the spins being detected. The proton gyromagnetic ratio Yh is about 4 times 
that of the carbon-13. Proton and carbon have different Ti and T 2 , which will also 
affect the SNR of the experiment. The maximum achievable SNR for the basic and 
signal enhanced J-editing methods are listed in Table 2. The J-editing methods with 
magnetization transfers use more pulses in the J-modulation buildup phase than the 
basic carbon detection methods. In inhomogeneous field, more pulses may result in 
larger errors and loss of signal due to imperfect pulses. The values listed in Table 2 
represent the maximum achievable sensitivities under ideal conditions. 

Table 2. Relative sensitivities of J-editing methods with/without magnetization transfers 

J-editing methods Relative sensitivity (SNR) 

Carbon detection J-editing ( Yc ) 5/2 =1 

Proton detection J-editing (y h ) 5/2 =32 

J-editing with l H -> 13 C magnetization transfer y h (y c ) 3/2 = 4 
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J-editing with 13 C-> l H magnetization transfer Y c (Yh) 3/2 = 8 

J-editing with *H ~* l3 C-± l H magnetization transfer (y h ) 5/2 = 32 



[0068] Two common techniques for magnetization transfer include insensitive nuclei 

enhanced by polarization transfer (INEPT) and Hartman-Hahn Cross Polarization. The 
INEPT technique transfers anti-phase proton spin magnetization to anti-phase carbon 
spin magnetization, or vice versa. See Morris, G.A. and Freeman, R., J. Am. Chem. 
Soc, 101, 760 (1979). This transfer is accomplished by applying a proton n/2 pulse and 
a carbon %/2 pulse. Using spin product operators, the magnetization transfer can be 
represented as: 

2I X S Z - (7c/2)I y ^ -2I 2 S Z — (tt/2)S x -> 2I z S y . 

[0069] Another common technique utilizes Hartmann-Hahn cross polarization (CP). See 

Hartmann, S R. and Hahn, EX., Phys. Rev:, 128, 2042 (1962). The proton and carbon 
spins are spin locked by the RF fields that satisfy the matching condition for equal 
precession rates in their respective rotating frames 

y h B* = y c Bj Eq.(6) 

where y h and y c are the respective gyromagnetic ratios for proton and carbon; Bj 1 and 
Bf are the respective RF field strength for proton and carbon. Magnetization transfer 
occurs at half of the J-coupling frequency. This can be used to transfer magnetization 
from proton to carbon, or from carbon to proton. 

[0070] To use the enhancement techniques in J-editing, it is necessary to design a pulse 

sequence having three parts: polarization, J-modulatiori, and detection. As noted above, 
the pulse sequences used in the polarization and J modulation periods may be referred 
to as a preparation pulse sequence. Embodiments of the invention may use various 
enhancement techniques (e.g., NOE, INEPT, CP) in the polarization and/or detection 
parts. The enhancement techniques involving NOE pulses have been described with 
reference to FIG. 8 and FIG. 9. 
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[0071] The following will describe a few exemplary pulse sequences that combine 

carbon-detected J-editing methods with the magnetization transfer techniques. 

[0072] FIG. 10 shows one embodiment of the invention using a pulse sequence for J- 

editing with 'H-^C magnetization transfer. As shown, the pulse sequence has three 
sections. In the first section ( ! H: (7c/2) x -x 0 -(7c)-T O -(7c/2)y, I3 C: (7c)^T 0 -(7t/2)), the proton 
magnetization is transformed into anti-phase proton magnetization through the J- 
coupling effect, and then is transferred into anti-phase carbon magnetization by the 
proton (7u/2)y pulse and the carbon 7c/2 pulse. The time duration x 0 is preferably equal to 
1/4 J- to maximize the magnetization transfer, where J is the average J-coupling constant 
of the sample. In the second section ( 1 H: x-(ti), 13 C: T-(7t)-T), the anti-phase carbon 
magnetization is refocused into carbon magnetization through the J-coupling effect. By 
shifting x, different amount of J-coupling effect on the carbon magnetization is 
achieved. The third section is the CPMG-like pulses. The carbon magnetization is 
detected in the form of a CPMG echo train. Multiple scans with different x values may 
be acquired for computing the relative carbon group abundances, as described above. 

[0073] FIG. 1 1 shows another embodiment of the invention involving a pulse sequence 

for J-editing with 13 C -> ! H magnetization transfer This pulse sequence also has three 
sections. In the first section ( 13 C (7t/2) x -T-(7i)-T, 1 H: t-(tc)), the carbon magnetization is 
transformed into anti-phase carbon magnetization through J-coupling. The amount of 
anti-phase carbon magnetization is changed when x is varied. In the second section 
( 13 C: (7t/2) y -Ti-(7r), l U: (7c/2)-ti-(7i)-ti), the anti-phase carbon magnetization is 
transferred into anti-phase proton magnetization and then refocused into proton 
magnetization. The time duration Xi is preferably a compromised value because carbon 
groups with different number of protons will behave differently. The third section of 
the pulse sequence is the CPMG-like pulses, during which proton spin echo signal is 
detected. The proton spin echo signal has different J-modulation effect when x is 
varied. 

[0074] FIG 12 shows a pulse sequence for J-editing with l H -» l3 C -> ! H magnetization 

transfer. This pulse sequence has four sections. In the first section ( l H: (7t/2) x -Xo-(7t)- 
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To-(7t/2) y , 1? C: (tc)-to-(7c/2) x ), the proton magnetization is transferred into anti-phase 
carbon magnetization. In the second section ('H: i-(7t), ,3 C: T-(tt)-T), the anti-phase 
carbon magnetization is modulated by a factor dependent on t. In the third section ^H: 
(7i/2)-ii-(7i)-ii, 13 C: (7r/2) x -Ti-(7r)), the anti-phase carbon magnetization is transferred 
into proton magnetization. In the fourth section, the proton magnetization is detected in 
the form of a CPMG echo train. 

[0075] In all examples of magnetization transfer described above, the n pulse after the 

time delay t may be replaced by a gated decoupling pulse as previously described in the 
gated-decoupling J-editing methods. The gated decoupling also destroys unwanted spin 
coherences. Therefore, similar pulse sequences using gated decoupling are expressly 
within the scope of the invention. 

Effects of Inhomogeneous Static and RF Fields 

[0076] In deriving Eqs. (1-3) above, the field inhomogeneity effects were neglected. 

However, when the pulses are inaccurate or the static field inhomogeneity is 
comparable to RF field strength, those equations are no longer correct, which is often 
the case in reality, especially in downhole environments. When the inhomogeneous 
field effects are taken into consideration, Eq. (1) becomes 

' S(T) = ^£ A if L VAV^ Eq. (7) 

i 

where j is the unitary imaginary number, V 0 is the total volume of the sample. The 
symbols L c 0 _, A c _ + , and A h 00 are matrix elements for rotation operations due to RF 
pulses. In spherical notation, the magnetization vector is defined as 

M + - M x +iM y 

M. -M x -iM y Eq. (8) 

M 0 =M Z . 

[0077] The rotation of the magnetization vector due to a RF pulse is represented by a 

matrix R, 
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where t p is the duration of the pulse. Matrix R is denoted by L for a tt/2 pulse and A for 
a 7i pulse. Furthermore, L 9 and A c denote matrices for carbon n/2 and 7i pulses, 

respectively; L h and A h denote matrices for proton k/2 and n pulses, respectively. The 

matrix elements L c 0 _ , A c _ + , and Ao, 0 in Eq. (7) are given by 



[1 - cos(Q c t° /2 )] + j sin(Q c t< /2 ) exp(j<|> c n/2 ) 



Eq(lO) 
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cos(Q h t*) . 



Eq.(ll) 



Eq. (12) 



[0078] The carbon Larmor frequency offset Aco^ is defined as Aa> c 0 = y c B 0 -o 0 ^, where y c 

is the gyromagnetic ratio of 13 C, B 0 the static magnetic field strength, and the RF 
angular frequency for carbon. The term ©J is defined as cof = y c Bj 12 where is the RF 
field strength for carbon. The carbon nutation frequency is , given by 
Q° - ^(Awo) 2 +(co|) 2 . The time durations of the carbon 7i/2 and 71 pulses are denoted by 
i c nl2 and t* , respectively, and the phases of those pulses are denoted by §\ l2 and ^, 
respectively. The proton related terms AcoJ , ©J 1 , Q h , t£ are defined similarly as those 
corresponding terms for oarbon. The integral in Eq. (7) integrates over the whole 
volume of the sample. 

[0079] The matrix elements L c 0 _ , A c _ + , and A h 0 0 can be theoretically computed after the 

static and RF fields are carefully mapped. Alternatively, these parameters can be found 
through calibration using some oil samples whose compositions are already known. In 
general, five parameters need to be found through calibration, which are jL° 0 ;_A c _ + dv, 
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jLVA-.+Ao^v , jLVA c . j+ (A h o, 0 ) 2 dv, j^A*! + (A h 0 , 0 ) 3 dv , and |L c 0 ^A%(A h 0>0 ) 4 dv . After 
these five parameters are determined, S(x) becomes a 4 th order polynomial of cos27n:J 1 . 

The abundance of different carbon groups can then be resolved by the means previously 
described for solving Eq. (1). 

[0080] Eq. (7) looks complicated. However, all complications come from the term A h 00 . 

If the proton n pulse is perfect everywhere in the sample, A h 0 0 becomes -1. The term 
[1(1 + A h 00 ) + X(l- A h 00 )cos27iJ i t] ni in Eq. (7) may then be simplified to [coslnhjp , which 
can then be taken out of the integral because it is not spatially dependent. With that 
simplification, Eq. (7) becomes Eq. (1) except for a constant factor |L C 0 _A C _ + dv that can 

be derived from calibration. 

[0081] Although it is difficult to have perfect proton k pulse everywhere in the sample in 

inhomogeneous fields, it is possible to avoid/minimize the J-modulation buildup after 
the imperfect proton % pulse, for example, using a gated-decoupling carbon-detection 
pulse sequence. The decoupling pulses should completely decouple the proton-carbon 
J-coupling such that Eq. (2) is correct even in inhomogeneous fields. Complete 
decoupling is achievable by using proper decoupling pulses or an RF transmitter 
coupled with a noise modulator to produce a broad band decoupling pulse. 

[0082] For carbon-flip proton-detection J-editing, Eq. (3) changes to the following form 

when the field inhomogeneity effects are taken into account: 

S(x) = if LV A h _, + (1 - A c 0 , Q )dv S ° £ A.H; (1 - cos 27TJ.T) . Eq.(13) ; 

Z^A 1 n 1 | 

i 

[0083] The term ^L h 0 _A h _ + (l- A c 00 )dv is a constant factor, which may be easily 

calibrated. The matrix elements L h 0 _, A h + , and A c 0 0 are analogous to L% _, A c _ + J and 
A h 0 0 that are given in Eqs.(10-12). 

[0084] FIG. 1 3 illustrate a method in accordance with one embodiment of the invention 

for obtaining J coupling information using J editing. As shown in step 1301, a series of 
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J-edited NMR measurements are acquired with different x delay times to the 
development of J modulations. These measurements may be performed with any of the 
pulse sequences illustrated in the present invention or with modified sequences that are 
analogous to those shown. A typical pulse sequence for this purpose may include two 
parts: the first part is a "preparation pulse sequence," which may include only a J-edit 
pulse sequence or additionally include a signal enhancement pulse sequence, and the 
second part is a "detection sequence." The J-edit pulse sequence permits J coupling to 
develop during a selected delay time x to modulate the detected signal amplitudes. The 
J-edit pulse sequence may include a pulse or a composite pulse to flip the coupled 
nuclear spins or a decoupling pulse to decouple the coupled nuclear spins. As noted 
above, the J-edit sequence may be used to develop J modulation for the first echo or for 
the first few echoes. The "detection sequence" may simply be a wait time (i.e., 
includes no pulse sequence) for the detector to acquire the first echo signal that resulted 
from the preparation pulse sequence. Alternatively, the detection sequence may include 
one or more n pulses to induce further echoes. When the detection sequence includes a 
plurality of n pulses, the overall sequence, including the preparation pulse sequence, is 
reminiscent of a CPMG pulse sequence. Such a CPMG-like sequence (i.e., a plurality 
of 7r pulses) may be used to collect a large number of spin echoes in a relative short 
period of time to improve the SNR. 

[0085] As shown in step 1302, the series of NMR measurements are then used to derive 

the J coupling information and/or the relative abundances of the carbon groups. This 
derivation may be performed by solving a series of linear equations that define the 
signal amplitudes as a function of the delay time x. Alternatively, this may be 
performed with a Fourier transformation, if the different x values are incremented by the 
same amount. The derivation involves correlating the signal amplitudes with the delay 
time x using a selected function that is based on the pulse sequence used in the data 
acquisition, e.g., Eq. (l)-(3), (7), and (13). 

[0086] The above described methods derive the J-coupling information to obtain the 

relative abundances of the carbon groups. In an alternative approach, the J-edit 
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spectrum (or the J-modulation curve as a function of delay time x) of a formation fluid 
may be deconvoluted into weighted contributions from a basis set of J-edit spectra (or J- 
modulation curve) recorded with standard samples having known compositions (e.g., 
methane, heavy oil, medium oil, light oil, and oil-based mud). The deconvolution may 
be accomplished by performing eigenanalysis based on the basis set of J-edit spectra or 
J-modulation curves. The relative weights of the basis set represent the relative 
abundances of the standard fluid types in the formation fluid. This alternative approach 
does not involve derivation of the J-coupling information. Instead, it directly provides 
relative abundances of various fluids in the formation fluid sample. This approach is 
not sensitive to magnetic field inhomogeneity and, therefore, no correction for the filed 
inhomogeneity is necessary. 

[0087] In some embodiments, the method may farther include a step for hydrocarbon 

typing based on the derived J coupling constants and the relative abundances of the 
carbon groups (Ai) (shown as step 1303). For example, if the J coupling constants 
include a population of 160 Hz, then the sample includes aromatics. Similarly, J 
, coupling of 125 Hz indicates the presence of aliphatic hydrocarbons. The relative 
contents of various hydrocarbons may be further classified based on the relative 
abundance Ai and nj. 

[0088] FIG. 14 shows a method in accordance with the invention for fluid typing. As 

shown in step 1401, a difference data set is obtained with proton detection. The 
difference is between two measurements, one with J modulation and the other without. 
The difference data set would include only those protons having 13 C nuclei attached to 
them. Proton signals due to protons not attached to 13 C nuclei are removed, including 
those signals from water. 

[0089] In step 1402 a hydrocarbon content (index) is derived from the difference data set. 

In deriving the hydrocarbon index, it is necessary to take into account that only about 
1 . 1% of carbons in nature are carbon- 13 species. 

[0090] In step 1403 the total hydrogen index is obtained with a conventional NMR 

measurement, which measures all signals from hydrogen species, be it from water or 



25 



PATENT APPLICATION 
ATTORNEY DOCKET NO. 20.2862 (R&O 09428. 106001) 

hydrocarbons. This total hydrogen index together with the hydrocarbon index can then 
be used to derive the relative abundance of water and hydrocarbons in the fluid sample. 
From the relative abundances of water and hydrocarbons, it is then possible to calculate 
the oil-water ratio or for fluid typing (shown as step 1404). 

[0091] Advantageously, the embodiments of the present invention provide convenient 

methods for determining relative abundance and types of various carbon groups, 
including quaternary C, aromatic CH, aliphatic CH, CH 2 , CH 3 , and CH 4 . Methods in 
accordance with the invention are based on J coupling and are less sensitive to magnetic 
field inhomogeneity compared to prior art methods. Because these methods are less 
sensitive to magnetic field inhomogeneity, embodiments of the invention may be used 
in a laboratory or in an environment where magnetic field homogeneity is difficult to 
obtain, for example, in a wellbore. 

[0092] The carbon group information obtained with methods of the invention may be 

used in many applications: measurement of hydrogen index of hydrocarbons in the 
presence of water, methane gas detection from the CH4 abundance; CH 2 /CH 3 ratio 
estimation, which is a good indicator of average chain length of the molecules; 
aromatics/aliphatics ratio determination; the abundance of aliphatic CH, which is 
related to branching and affects wax formation; and quaternary C estimation as an 
indication of asphaltene, to name a few In addition, carbon group abundances thus 
determined may provide useful inputs for molecular weight analysis based on diffusion 
or T 2 distributions. - 

[0093] While the invention has been described with respect to a limited number of 

embodiments, those skilled in the art, having benefit of this disclosure, will appreciate 
that other embodiments can be devised which do not depart from the scope of the 
invention as disclosed herein. Accordingly, the scope of the invention should be limited 
only by the attached claims. 
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